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Abstract 

It is shown that contrary to a recent claim in the literature0, the domain walls made of 
scalar field are diamagnetic due to presence of massless fermionic modes on the wall. The 
diamagnetism vanishes at high temperature. Thus the domain walls could produce no effect 
on a primordial magnetic field in the early universe. 



In a recent paper IwazakiQ has considered the magnetic properties of the gas of massless 
fermions bounded to a domain wall in a theory with spontaneously broken discrete symmetry. 
The result is claimed to be that the free energy of the gas in magnetic field with the strength B 
contains a negative term proportional to —5^/^, thus that the gas should exhibit spontaneous 
magnetization. If this were true, the domain walls in the early universe could be a source of 
the primordial magnetic field. 

However, the calculations in Ref. 1 are not quite correct, and the purpose of this note 
is to present a correct calculation of the magnetic properties of the gas of massless fermions 
on the wall at zero temperature as well as at finite temperature. The result of the present 
calculation is that in fact the term in the free energy proportional to B^^"^ is positive and 
goes to zero at high temperature ^ eB. Thus the domain walls are diamagnetic, and 
they can not generate a primordial magnetic field in the early universe. 

It can be reminded that a domain wall, located for definiteness perpendicular to the z axis, 
corresponds to transition of a scalar field (f){z) between two minima vi and V2 of the scalar 
field potential: (f){z = — oo) = fi, (j){z = +00) = V2- Furthermore, if there is a fermion field 
with a Yukawa coupling to such that the fermion mass term m(0) changes sign between vi 
and f 2 such fermion has a zero mode with respect to motion in the z directionB . The fermion 
wave function corresponding to this zero mode thus splits into the product ip{t,x,y, z) = 
x(t,x,y) u{z) with the scalar function u{z) given by u{z) = const ■ exp (± J^mdz) and 
the spinor x being an eigenvector of the Dirac matrix ^73: i'jsx = where the upper 
(lower) sign corresponds to 'm{v2) negative (positive). The spinor condition on x leaves only 
two components independent, thus the motion of the zero-mode fermions along the wall is 
described by the (2+1) dimensional Dirac equation for a two-component massless spinor. 
Clearly, the solutions to that equation (on-shell particles) have one degree of freedom for 
fermions and one for antifermions. In the absence of an external magnetic field the spectrum 
of the free Dirac operator is obviously labeled by the two components of the two-dimensional 
momentum k: E{k) = ±k (with k = |k|). In the presence of a uniform magnetic field B 
orthogonal to the wall, the spectrum of the Dirac operator with the electric charge e is 
foundS following the textbook method of LandauS and is parametrized by one integer 
quantum number n > with the degeneracy eB/{2Tx) per unit area of the wall and the 
energy given by En = ±V2 e B n. 

The existence of the non-zero modes of the motion in the z direction, either non-localized 
or localizedS at the wall, can be ignored in calculating the magnetic and thermodynamic 
properties of the gas of the zero-mode fermions as long as the characteristic energy scale in 
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the problem, i.e. ^eB and/or the temperature T, is much smaller than the mass gap given 
by the mass of the fermions in a spatially uniform vacuum state v\ or v^- Throughout this 
letter this is assumed to be the case and the contribution of non-zero modes to the partition 
function is completely ignored. 

The standard construction of the physical states of the fermion field starts with defining 
the vacuum state, where the negative energy states are occupied and the positive energy 
ones are vacant. In the considered here system of the two-dimensional massless fermions 
this leaves a (^2)^ degeneracy of the vacuum states corresponding to the two possibilities 
for the occupation number at n = and the total number of such states is = x Area. 
The non-vacuum states, as usual, correspond to fermions: filled states with positive energy, 
and to antifermions: holes in the states with negative energy. Accordingly, the free energy 
per unit area of the fermionic gas at a temperature T = in the magnetic field B can be 
written in the following form 

F = F_\F^\F^\E,^ , (1) 

where 

F± = -/?-^ ^ ^ In (1 + e-/3v^^) (2) 

2^ n=l 

is the free energy associated with the real gas of fermions and antifermionsQ. The term Fq 
in eq.(P is associated with the degeneracy of the vacuum state: 

Fo = -r ' If 1^ 2 . (3) 

Finally, i?vac is the energy per unit area of the wall of any of the degenerate vacuum states^. 
Clearly, i^vac gives the free energy at T = 0. 

Ignoring the divergence of the sum, one might write the vacuum energy as the sum over 
the energies of the occupied negative-energy eigenstates: 

^vac = -i^ E v/275;i=-—v/275C -- - 0.04679 (e 5)^/' , (4) 



^The fermion number can be freely exchanged with the non-zero modes, hence there is no chemical 
potential in eq.(^0. 

^It is this vacuum energy term which is missing in the corresponding calculation of Ref. 1. Also the 
Euler-Maclaurin expansion for the sum in eq.(H) is unjustifiably truncated, which has lead to a wrong result 
there. In fact this expansion is of little help in calculating the sum because of the root singularity of the 
summand at n = 0. 
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where the finite answer for the divergent sum is written in terms of the standard ({s) function 
analytically continued below its pole at s = 1. It turns out that this formal manipulation 
gives the correct expression for the dependence of the vacuum energy on the field strength B, 
while the overall vacuum energy is of course infinite. In view of this divergence a somewhat 
more elaborate consideration of the vacuum energy is due, and this consideration, presented 
in the following lines also introduces an appropriate method for calculating the sums involved 
in this problem^. 

First, in order to make tractable the sum over the infinite number of negative energy 
states, it has to be regularized. A gauge invariant regulator factor should depend on a gauge 
invariant quantity, which in this case is naturally the level energy itself. We chose here the 
regulator factor in the exponential form: exp{—eE^) with e being the regulator parameter. 
Thus we write the regularized vacuum energy in the form 

eB °° 

^i:^ = -^ E v/275;^exp(-e2ei?n) . (5) 

n=l 

The method of summation to be used here is based on the identity valid for any smooth 
function f{x): 

/-oo poo 

E /(^) = / /(^) E 5{x — n) dx = E / f (x) exp{2 IT i m x) dx , (6) 

n=l ^ ?i=— oo m=— oo ^ 

where S is any number such that < S < 1. Notice that the sum in the second expression 
goes over n from — oo to +oo. The identity is still valid since the terms with n < are 
identically zero. Applying this identity to the sum in eq.@ and noting that in this case in 
fact S can be set equal to zero, since the term in the sum with n = is vanishing, one can 
write 

-^vac = E / V2e B X exp (— e 2e B x) exp(2 vri ma;) dx . (7) 

2ti ^_ _ JO 



For each m the integral in this equation is clearly that for the F function: F(3/2). The term 
with 171 = diverges in the limit e — > 0, while all the rest terms of the sum are finite. Thus 
writing separately the singular term and grouping together the terms with symmetric values 
of m one finds for e ^ 

Svre'*/^ 27r An 



■^I am thankful to Igor Aleiner for reminding me this summation technique, which in fact was used by 
Landau in similar problems and which originally goes back to Poisson. 
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The singular term in this expression does not depend on B and is exactly equal to the energy 
of the free vacuum regularized in the same way: 



411(0) (9) 

while the second term in eq.(|D, which is finite and describes the dependence of i?vac on the 
magnetic field is exactly equal to that in eq.(|D due to the identity ({—1/2) = — C(3/2)/(47r), 
which is a consequence of a general relation0 between ({s) and ({1 — s). 

The term with S^/^ describing the dependence of the vacuum energy on B in eq.(^ or 
in eq.(|^) is positive. This corresponds to diamagnetism of the wall at zero temperature with 
a singular at B ^ diamagnetic susceptibility. The temperature dependence is given by 
the terms F± and Fq in the free energy in eq.(|l|). Following Ref.l, we consider here the free 
energy at a large temperature, i.e. in the limit (3 = (3 \/l e B <^ 1 (but still the temperature 
is much less than the mass gap in order to ensure the irrelevance of the non-zero modes). 
According to eqs.(0 - the temperature dependent part of the free energy is given by 

.1 eB 

Ln=l 



ln(l + e-'^^) + - ln2 



(10) 



In order to find the expansion of the sum in this expression in powers of (3 at small (3 we 
apply to the sum the identity (|^). Writing separately the Fourier harmonic with m = and 
grouping together the harmonics with symmetric m, one finds 



Eln(l + 

n=l 

POO - '"'^ POO 

/ ln(l + e-^^) ^ / ln(l + e-^v^) (e^'^^'"^ + 6-2'^^'"^) rfx . (11) 

The first integral in the limit 5 — +0 gives the singular in /3 term | which upon 
substitution in eq.(0) reproduces the field independent free energy of a free fermion gas. 
The rest of the terms of the expansion in powers of /3 are obtained by the Taylor expansion 
of the integrand in the second term in eq. ([Tl|) with subsequent separate integration and 
summation over m in each term of the expansion in (3. The parameter 5 can be set to zero 
before integration in all the terms of the Taylor expansion except the first one, where keeping 
5 small but finite ensures the convergence of the sum over m. The calculation of this first, 
constant, term is thus reduced to the formula 

lim y / fe2"^'"^' + e-2"^™^) c/a; = — (12) 
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and the higher terms of the Taylor expansion in (3 are found from the generic formula 
£ [ x^l^ (e—^ + e-—^) dx = (2.)-/^ ^ (| + l) C (| + l) (13) 



m=l 



for p > 0. (Notice that eq. can also be found by taking in this formula the limit p — » 0.) In 
deriving these formulas the oscillating behavior of the integrand is dealt with by temporarily 
introducing the damping factor exp(— ex) and taking e — > +0 in the end: this is a legitimate 
procedure, since it does not alter the original integral in eg. (pJ]) . 

It is satisfying to see, using eq.(|12[), that upon substitution in eq.(p!OD the constant, [3 
independent, term in the sum of eq.(|TO[) exactly cancels the contribution Fq to the free energy 
of the vacuum degeneracy factor. Thus, as expected on general grounds, no term linear in 



B arises in the free energy at any temperature. Therefore, using the formula in eq. ([T3|) for 
p = 1 one finds in the high temperature limit the expansion for the thermal part of the free 
energy up to the first B dependent term as 



P^^F_^F, = -l^-f-^:^^^-m^O{,H.Bf) . (14) 



One can readily see that the term with B^/"^ in this expression exactly cancels in the sum 
in eq.(|^) the corresponding contribution of the vacuum energy, given by eq.(^. Thus one 
concludes that the singular at 5 ^ diamagnetic behavior of the fermion gas at the wall 
vanishes at high temperature, as could be expected on general grounds. 

It is a simple exercise to calculate numerically the sum in eq . ([T0|) for arbitrary /3. However 
in view of the absence of any dramatic phenomena, like the previously claimedfl ferromag- 
netism of the fermion gas, this calculation is of a little practical interest. 

As discussed above, the calculation here is restricted to the situation where the tem- 
perature is much smaller than the mass gap for the fermions. If this restriction is relaxed, 
one has to solve the problem for the full spectrumS of the fermion scattering states in the 
presence of the domain wall in (3+1) dimensions. In this case however, the effect of the wall 
on the magnetic properties of the field system is small, and also no dramatic effects can be 
expected. 

Thus we conclude that the domain walls, had they existed in the early universe in spite 
of their undesirable cosmological effects pointed out in the literature©, would not have 
produced any significant impact on a large-scale primordial magnetic field. 

I am thankful to I. Aleiner, L. Glazman, A. Larkin and A. Vainshtein for enlightening 
discussions. This work is supported, in part, by the DOE grant DE-AC02-83ER40105. 
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